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ABSTRACT 


A  hybrid  finite  element  cum  T-matrix  technique  is  introduced  for  the 
analysis  of  problems  involving  acoustic  wave  scattering  by  shells  of  revo¬ 
lution  immersed  in  water.  The  concept  of  a  mechanical  impedance  or  recep- 
tance  is  used  from  the  theory  of  vibrations  to  provide  a  relationship 
between  the  velocity  of  the  shell  surface  and  the  pressure  acting  upon  it. 
This  relationship  is  introduced  in  the  Helmholtz  integral  relations  for  the 
incident  and  scattered  fields  and  a  T-matrix  for  the  scattering  problem  is 
derived.  The  mechanical  impedance  is  calculated  by  a  conventional  FEM 
technique.  Numerical  results  are  obtained  for  spherical  shells  as  a  check 
on  the  numerical  procedure  and  also  for  a  finite  capped  cylindrical  shell 
for  waves  incident  along  the  rotational  axes  of  symmetry.  The  comparison 
with  the  experimental  results  of  Dr.  S.Numrichand  Dr.  L.  Dragonette  of 
the  Naval  Research  Laboratory  is  excellent. 
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INTRODUCTION 

1  2 

In  chis  report  a  hybrid  method  employing  the  T-matrix  approach  * 
along  with  a  finite  element  formulation  is  introduced  for  the  analysis 
of  problems  involving  the  scattering  of  sound  waves  by  elastic  shells 
immersed  in  water.  It  is  shown  that  by  using  the  concept  of  a  mechanical 
impedance,  we  can  avoid  the  use  of  integral  representations  in  the  Interior 
of  the  scatterer  or  the  use  of  fluid  elements  in  the  exterior  and  3-D 
elements  in  the  interior  of  the  shell.  In  other  words,  this  method  ex¬ 
ploits  the  best  features  of  both  the  T-matrix  method  and  the  FEM.  The 
latter  of  course  can  successfully  handle  several  geometries  of  practical 
interest . 

The  impedance  of  the  shell  is  calculated  and  it  provides  the  necessary 
relationship  between  the  velocity  on  the  shell  surface  and  the  pressure 
acting  on  it.  Thus  the  number  of  unknowns  in  the  scattering  problem  is 

j 

reduced  by  one.  The  integral  representation  of  the  scattered  field  and 
the  extincting  theorem  are  then  used  in  the  conventional  manner  to  yield 
a  T-matrix  for  the  problem.  Another  byproduct  of  the  receptance  calcula¬ 
tion  is  that  the  free  vibration  frequencies  of  the  shell  are  also  generated 
many  of  which  can  be  identified  with  the  maxima  and  minima  of  the  spectrum 
of  the  scattered  pressure  field. 

Previous  calculations  for  acoustic  wave  scattering  from  elastic  shells 
of  revolution^  experienced  numerical  difficulties  as  the  shell  thickness 
decreased  and  the  frequency  increased.  Moreover  the  formulation  used  • 

invoked  integral  representations  in  the  interior  of  the  shell  region,  this 
resulted  in  very  complicated  expressions  for  the  T-matrix  which  depended 
on  5  or  6  different  matrices.  Some  of  these  matrices  are  quite  ill  condi¬ 
tioned.  With  the  present  scheme,  the  scattering  problem  is  a  purely  scalar 
one  and  hence  much  better  conditioned  numerically. 
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FORMULATION  OF  THE  SCATTERING  PROBLEM 

Consider  a  shell  of  revolution  bounded  by  the  surface  S  with  a 
piecewise  continuous  outward  normal  n  that  is  immersed  in  water.  The 
acoustic  field  in  water  is  described  by  the  scalar  potential  $.  We  con¬ 
sider  a  time  harmonic  plane  wave  of  frequency  tu  incident  on  the  shell. 
Without  loss  of  generality,  the  plane  of  incidence  can  be  considered  to 
be  the  x-z  plane.  Further  the  time  dependence  of  all  fields  is  exp(-lut) 
and  this  is  suppressed  for  brevity.  The  total  field  outside  the  shell 
is  given  by 

<|>(r)  -  $Q(r)  +  $g(r)  ;  r  e  V  (1) 

where  4>  and  $  are  the  incident  and  scattered  fields  respectively, 
o  s 

The  free  space  Green's  function  in  the  fluid  is  given  by 

g(r,r')  -  expikf |r-r' |/ |r-r’ |  (2) 

where  k^  *  w/c^  is  the  wavenumber  in  the  fluid,  cf  is  the  sound  speed 
in  the  fluid  and  r  and  r'  are  respectively  the  field  and  source  points. 

The  scalar  potential  satisfies  a  homogeneous  wave  equation  in  V  and 
so  also  does  g  except  at  r  -  r ' .  The  integral  representation  for  such 
a  system  is  known  to  be 


JL 

4ir 


/  u;n'-v'g(?,?') 
S 


g(r  ,r ')n'  }dS' 


-*0(r)  ;  r  t  V 

;  *  e  v 


(3) 


The  first  form  of  Eq.  (3)  for  r  t  V,  is  known  as  the  extinction  theorem 
and  the  second  form  of  Eq.  (3)  for  r  e  V  is  the  integral  representation 
of  the  scattered  field.  The  surface  quantity  $+  -  (iuJOf)'1?^..  where 


P+  is  the  surface  pressure  and  p^  is  the  mass  density  of  the  fluid.  The 
surface  quantity  n*?+$  •  n*V^_  where  V+  is  the  fluid  particle  velocity  on  S. 
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(9) 
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1 r  /  (M-VRe,('4mo  "  Re*W**V}dS 


The  surface  pressure  and  normal  component  of  the  particle  velocity  are 
related  through  the  mechanical  receptance  of  the  shell  as  follows. 

We  expand  the  normal  component  of  the  velocity  on  the  shell 


surface 


"  *  V  -  l  Cn  ‘Jn 


where  qn  are  any  set  of  basis  functions  associated  with  the  shell.  In 
the  next  section,  we  shall  identify  qR  with  the  eigenvectors  of  the  shell. 
The  pressure  on  the  surface  of  the  shell  is  expanded  as  shown  below 

♦+  -  E  B„  <12> 

n 

The  mechanical  receptance  of  the  shell  is  the  ratio  of  the  surface  ve¬ 
locity  to  the  surface  pressure  and  In  the  next  section  we  show  that 

C  =•  I  R  ,  B  ,  (13) 

n  L,  nn  n' 

n 

where  R  is  the  receptance  matrix.  Substituting  Eqs.  (11)-(13)  in  (9)  and 
(10)  we  have 

a  *  -i  l  d4) 

n  i  iui  ti  *n  n 

n 


fn  •  1 


where  the  super  index  n  -*•  (i,m,a)  and 


<W(0u,lWV>  *  Z7  l„  .[  {("'7Re  W 

n  s 


!  Vn"  ‘  Rn,n,,Ee*««o  }qn"dS 


(16) 
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Equations  (14)  and  (15)  can  now  be  formally  solved  to  yield 
f 


J,  T  *  i 
,  nn  n 


where 


T  -  -  Q(Reii»,q)  [Q(Ou^,q)  ] 


-1 


(17) 


(18) 


is  the  T-matrix  of  the  scatterer. 

The  following  properties  of  the  T-matrix  resulting  from  reciprocity 
and  the  conservation  of  energy  can  be  shown 

T  *  TC  (symmetry)  (19) 

TT*  -  -  Real(T)  (20) 

The  above  properties  can  be  used  as  a  check  of  the  numerical  computations 
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CALCULATION  OF  MECHANICAL  RECEPTANCE 

To  calculate  the  receptance  the  shell  theory  of  elasticity  is  used 
in  conjunction  with  a  finite  element  approximation.  Consider  a  shell  of 
revolution  of  thickness  2h.  The  variable  s  is  a  length  coordinate 
measured  on  the  mid-surface  oi  the  shell  (see  Fig.  1).  Let  u  and  w  be  the 

normal  and  tangential  displacements  of  the  shell  surface  and  8  is  a  meridi- 

2  3 

onal  deflection.  According  to  the  assumptions  of  shell  theory  *  ,  the 
strains  and  displacements  are  related  as  follows 

3u  w 
3s  R 

(wcos$  +  usin$)/r 

d8  d  .  ,  , 

-  dl  +  di(u/R) 

-  sin<j>  (8-u/R)  /r 
dw 

-  8 
ds 

where  R  is  the  principal  radius  of  curvature  in  the  u-w  plane  and  r  and 
^  are  measured  as  shown  in  Fig.  1.  In  Eq.  (21),  eg  and  are  membrane 
strains,  xg  and  are  bending  strains  and  Y  is  a  transverse  shear  strain. 
The  strain  energy  of  the  shell  can  be  written  as 


It  =*  y  /  eT  D  e  dS  -  /  fT  u  dS  (22) 

where 

(23) 

is  the  displacement  vector  and  f  is  the  vector  of  the  external  force  and 
0  is  the  elastic  matrix 


0 


(24) 
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In  Eq.  (25),  E  is  the  Young's  modulus  and  v  is  the  Poisson's  ratio. 

4 

Through  the  usual  finite  element  approximation  process  ,  the  dis¬ 
placements  at  any  point  on  the  surface  of  the  shell  are  written  as 
D(r)  -  N±U1  +  K1+l  Ui+1  ;  r  e  Et  (26) 

where  are  shape  functions  and  i  refers  to  the  i-th  node  and  is  the 
i-th  element.  and  U^+^  are  the  nodal  values  at  the  edges  of  the  element 

E^.  We  note  that  for  axisynmetric  problems,  there  is  no  dependence  of  the 
fields  on  the  azimuthal  angle. 

Using  Eq.  (26)  in  (21)  and  (22),  the  element  stiffness  matrix  can 
be  written  as 


r+1  T 

K  -  /  B^DBi+]dS 


where 


„  1-n  .  „  1+n 

Ni  *  ~r and  Mi+i  *  ~r 


where  ' n'  is  measured  along  the  curved  length  of  the  element. 
The  matrix  B  is  given  by 
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The  element  mass  matrix  can  be  derived  from  the  expression  for  the 
kinetic  energy  which  after  Integration  in  the  thickness  direction  and 
the  azimuthal  direction  can  be  written  as 

2 

T  -  y(ph)  /  [u2  +  w2  +  yj  02<J>  ]2rrrds  (29) 

where  p  is  the  mass  density  of  the  shell.  The  last  term  in  the  integrand 
of  Eq.  (29)  is  the  rotary  inertia  term. 

The  minimization  of  the  Lagrangian  leads  to  the  following  matrix 
equat ion 

K{Q>  +  M{Q}  -  0  (30) 

where  (Q>  is  the  nodal  vector  of  modal  values  of  the  displacement  and  M 
is  the  mass  matrix  that  results  when  Eq.  (26)  is  substituted  in  (29)  and 
the  Lagrangian  is  minimized.  If  one  assumes  that  the  time  dependence  of 
{Q>  is  of  the  form 

iw  t 

{Q}  ■  qne  (31) 

then 

K{qn)  »  -a>n2M{qn}  (32) 

where  is  the  n-th  eigenvalue  and  tq^}  is  the  corresponding  mode  shape. 
Since  the  eigenvectors  are  arbitrary  up  to  a  multiplicative  factor,  it 
is  convenient  to  make  use  of  the  fact  that  {q^}  are  orthogonal  with 
respect  to  the  mass  matrix  and  normalize  the  {q^}  such  that 

{q  }T  M { q  .  }  -  <5  ,  (33) 

n  ’n  nn 

Once  the  eigenvalues  and  eigenvectors  are  determined  for  the  free 
vibration  of  the  shell,  now  the  forced  vibration  problem  is  considered. 

The  boundary  condition  on  the  shell-fluid  interface  is  that  the  pressure 
and  the  normal  component  of  the  velocity  are  continuous.  The  surface 
value  of  the  pressure  is  related  to  and  the  equivalent  nodal  forces 
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can  be  written  as 


{F}  -  (f}e"iU)t  -  l  Bn{f 

n  ti 


(34) 


where  the  {fR}  are  derived  from  Eq.  (12).  The  vibration  problem  now  is 


K(Q}  +  M(Q}  -  (F) 


(35) 


The  nodal  displacements  are  expanded  in  terms  of  the  eigenmodes  as 

{Q}  -  l  C  (q  } 
u  n  nn 


(36) 


Substituting  Eqs.  (34)  and  (36)  in  (35),  we  obtain 


nn 


2  2 

W 

n 


(3 


This  is  the  form  of  R  used  in  the  definition  of  the  Q-matrix  in  Eq.  (16) 
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RESULTS  AND  DISCUSSION 

The  first  chock  of  this  method  was  done  by  performing  calculations  for 
a  spherical  shell.  The  top  hemisphere  of  the  shell  was  successively  divided 
into  10  and  then  20  elements.  Both  symmetric  and  anti-symmetric  boundary 
conditions  were  applied  at  the  equator  and  the  eigenvalues  and  eigenvectors 
were  generated  by  the  usual  finite  element  procedure.  It  was  judged  that 
20  elements  wore  sufficient  to  assure  convergence.  The  eigenmodes  of  the 
sphere  fall  into  two  broad  classes,  the  so-called  bending  and  membrane  nodes. 
In  Table  I,  the  natural  frequencies  for  the  spherical  shell  computed  with 
10  and  20  elements  are  compared  with  analytical  results. 

For  the  spherical  capped  cylinder  where  length  (22.)  to  diameter  (2a) 
ratio  is  1.6,  (see  Fig.  '*)  ,  35  elements  were  used  for  the  upper  half  of  the 
shell.  In  Table  II  the  natural  frequencies  of  the  capped  cylinder  are  tabu- 

GJ  i2 , 

lated.  For  the  scattering  problem  a  frequency  range  0<-= — <15.0  was  considered 

2Cf 

and  44  eigenmodes  were  used  in  the  calculation  of  the  receptance.  Out  of 
these,  25  can  be  identified  to  be  predominantly  bending  modes.  The  symmetry 
of  the  T-matrix  was  better  than  1/2%  in  the  entire  frequency  range 
0  <  k^i/2  <  15.0.  Only  the  case  of  waves  incident  along  the  rotational 
axes  of  symmetry  was  considered.  The  frequency  response  of  the  scattered 
pressure  field  at  distances  far  from  the  shell  are  plotted  in  Figs.  3,  4 
and  5  for  3  different  angles  of  observation.  In  Fig.  3,  a  portion  of  the 
calculated  back  scattered  spectrum  is  compared  with  the  experimental 
results  provided  to  us  by  Dr.  S.  Numrich  and  Dr.  L.  Dragonette  of  the 
Naval  Research  Laboratory,  Washington,  D.C.  The  agreement  is  excellent. 

The  experiment,  which  was  of  Che  pulse-echo  type  had  a  sampling  interval 
of  Ak^a  =  0.1  and  hence  the  sharp  resonance  at  k^a  =  9.3  which  is  seen 
in  the  calculation  was  missed  while  performing  an  FFT  from  the  time  to  the 
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frequency  domain.  The  experiment  is  very  well  calibrated  and  the  maxi¬ 
mum  difference  in  the  amplitude  is  less  than  3db. 

In  conclusion,  we  have  demonstrated  a  practical  technique  for  the 
study  of  acoustic  wave  scattering  by  thin  elastic  shells  of  revolution. 

Due  to  the  use  of  the  finite  element  method  for  the  calculation  of  the 
shell  receptance,  practical  shapes  and  features  such  as  stiffeners  can 
be  taken  into  consideration.  The  results  from  this  new  technique  are 
far  superior  to  that  resulting  from  the  use  of  a  full  elasticity  solution 
for  the  shell.  The  approach  is  better  conditioned  numerically  and  is  ex¬ 
pected  to  yield  stable  results  for  aspect  ratios  greater  than  that  by 
Che  conventional  solution.  It  now  remains  to  extend  the  calculation  to 
oblique  angles  of  incidence  with  respect  to  the  rotational  axes  of  symmetry. 

I 

1 
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TABLE  I 


Natural  Frequencies  ft  of  Spherical  Shell 
in  Nondimens ional  Units  ft  » 


10  elements 

20  elements 

Analytical 

b! 

0.2905  x  10~3 

0.6774  x  10~4 

0.0 

b2 

0.7221 

0.7245 

0.7259 

b3 

0.8756 

0.8814 

0.8855 

b4 

0.9786 

0.9858 

0.9939 

b5 

1.107 

1.108 

1.212 

b6 

1.295 

1.274 

1.296 

b7 

1.561 

1.495 

1.528 

b8 

1.915 

1.771 

1.818 

b9 

2.364 

2.100 

2.167 

b10 

2.914 

2.479 

2.570 

bll 

3.574 

2.904 

3.026 

b12 

4.351 

3.374 

3.532 

ao 

1.78713 

1.78713 

1.78713 

al 

2.187 

2.188 

2.189 

a2 

2.990 

2.988 

2.989 

a3 

3.981 

3.966 

3.963 

a4 

5.056 

5.004 

4.993 

a5 

6.158 

6.070 

6.044 

a6 

7.312 

7.154 

7.105 

a7 

8.510 

8.253 

8.172 

a8 

9.757 

9.367 

9.242 

a„ 

10.50 

10.31 

9 

a10 

11.64 

11.39 

all 

12.80 

12.46 

a.  „ 

13.98 

13.54 

12 

TABLE  II 


Natural  Frequencies  0  of  Finite  Capped  Cylinder , 
4/2a-1.6  in  Non-dimensional  Units  (fi  •  /i.R2<i>2/E) 

Bending  Modes  Membrane  modes 


bl  s 

0.0 

ao 

5.387 

b2  ” 

2.639 

al 

7.522 

b3  - 

3.213 

a2 

8.904 

b4  " 

3.696 

a3 

11.45 

b5  * 

3.854 

a4 

14.33 

- 

4.038 

a5 

17.08 

b7  * 

4.331 

a6 

20.13 

b8  * 

4.731 

a7 

23.11 

b9  ” 

5.207 

a8 

26.09 

b10  “ 

5.893 

a9 

32.23 

bll  “ 

6.53 

a10 

35.31 

b12  ' 

7.343 

all 

41.55 

b13  ' 

8.257 

a12 

44.7 

b14  " 

9.243 

b15  * 

10.29 

b„  ,  • 

11.48 

16 

b17  - 

12.71 

b18  “ 

14.01 

b,  _  * 

15.4 

19 

b__  * 

16.84 

20 

b21  * 

18.36 

b22  “ 

19.8 

b23  ’ 

21.55 

b24  " 

22.24 

FINITE  CYLINDRICAL  SHELL  IN  WATER 
FORM  FUNCTION  vs  FREQUENCY 


Fig.  3 


Backscattered  form  function  as  a  function  of  non- 
dimensional  wavenumber  for  normal  incidence  (9  >0, 
9*180°) . 


non-dimensionai  ».ve„u,ll,r 


